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Abstract. We study the effect of vaccination on robustness of networks against 
propagating attacks that obey the susceptible-infected-removed model. By extending 
the generating function formalism developed by Newman (2005), we analytically 
determine the robustness of networks that depends on the vaccination parameters. 
We consider the random defense where nodes are vaccinated randomly and the degree- 
based defense where hubs are preferentially vaccinated. We show that when vaccines 
are inefficient, the random graph is more robust against propagating attacks than 
the scale-free network. When vaccines are relatively efficient, the scale-free network 
with the degree-based defense is more robust than the random graph with the random 
defense and the scale- free network with the random defense. 
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1. Introduction 

Many real networks, such as the WWW, the Internet, social and biological networks, 
have complex connectivity. A property shared by many types of networks is the 
scale-free (SF) degree distribution pk oc k~'^, where k is degree (number of links 
connected to a node) and pk is the fraction of nodes with degree k [1]. Functions 
of networks crucially depend on the degree distribution and other structural properties 
of networks [21 El 13 El El El E] • 

An important property of networks is robustness against the removal of nodes 
caused by failures or intentional attacks. Albert et al. studied the robustness of networks 
against two types of attacks: random failure, where nodes are sequentially removed with 
equal probability, and intentional attack, where hubs (i.e., nodes with large degrees) are 
preferentially removed [9]. They showed that SF networks having 7 < 3 are highly 
robust against the random failure in the sense that almost all nodes have to be removed 
to disintegrate a SF network. However, SF networks are fragile to the intentional attack 
in the sense that the network is destroyed if a small fraction of hubs are removed. 
These results have been analytically established [lOl [HI [12] . The robustness of networks 
against other percolation-like processes such as the betweenness-based attack [13] and 
degree-weighted attacks [H] have also been studied. 

In fact, attacks to networks may occur as propagating processes. Examples 
include computer viruses [151 (HI (13 • The (first) critical infection rate above which 
a global outbreak occurs has been obtained for propagating processes such as the 
susceptible-infected-removed (SIR) and susceptible-infected-susceptible (SIS) models. 
In particular, for SF networks with 7 < 3, any positive infection rate induces a global 
outbreak [21 [3l [H [3, [6l [3, [HI [13, [HI [19] (also see a recent paper [20] in which it is shown 
that the critical infection rate of the SIS model on SF networks vanishes regardless of the 
value of 7). Another (second) critical infection rate, above which the network remaining 
after all the infected nodes are deleted (we call it residual network) is disintegrated, 
also exists and is larger than the first critical infection rate. The residual network 
is important because a second epidemic spread may occur on it [211 [221 [231 [21] • By 
using the generating functions, Newman derived the second critical infection rate for 
uncorrelated networks to show that the second critical infection rate is positive even for 
7 < 3 [21]. In particular, the global outbreak of a second epidemic requires a larger 
infection rate than that of a first epidemic [21] . 

Defense (i.e., immunization) strategies for networks to contain epidemics are also 
of practical importance. Epidemics can be efficiently suppressed by various defense 
strategies such as target immunization [23 [26] , acquaintance immunization [26l [23, [23 
[29] , and graph partitioning immunization [30]. Such a defense strategy raises the first 
critical infection rate. In contrast, we study the effectiveness of a defense strategy in 
enhancing the robustness of networks, i.e., increasing the second critical infection rate. 
Using the generating functions [211 [31], we formulate the effects of defense strategies 
on the robustness of networks. We determine the two critical infection rates for the 
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Figure 1. Propagating attack on networks (a) without and (b) with a defense strategy. 
T and T/v are the probabihties that an infected (and eventuahy removed) node infects 
an unvaccinatcd and vaccinated neighbor in a propagating attack, respectively. 

following three combinations of networks and defense strategies: random graph with 
the random defense, in which nodes are vaccinated in a random order, uncorrelated 
SF network with the random defense, and uncorrelated SF network with the degree- 
based defense, in which nodes are vaccinated in the descending order of the degree. 
When vaccines are inefficient, the random graph is more robust than the SF network in 
the sense that the second critical infection rate is larger. When vaccines are relatively 
efficient, the SF network with the degree-based defense is more robust than the other 
two cases. We also discuss the optimization of the vaccine allocation under a trade-off 
between the number and efficiency of vaccines available. 

2. Model 

We adopt the SIR model, which is a continuous-time Markov process on a given static 
network with nodes, as the propagating attack. Each node takes one of the three 
states: susceptible, infected, and removed. All nodes except a randomly selected seed 
node are initially susceptible. The seed node is initially infected. Susceptible nodes get 
infected at a rate proportional to the number of infected neighbors. In other words, if 
a susceptible node is adjacent to an infected node, the susceptible node gets infected 
with probability f3At within short time At. An infected node becomes removed at a 
unit rate, i.e., with probability At within short time At, irrespective of the neighbors' 
states. 

Following the methods employed in previous literature [311 [32], we analyze the 
bond percolation to estimate the transition points, the final size of global outbreaks, 
and that of residual networks of the SIR model. We assume that the disease is 
transmitted from an infected node to a susceptible node with probability T = T{P), 
where T = /[I — exp(— /3r)]P(r)dr, r is the infectious period during which the infected 
node is infected, and -P(t) is the distribution of r [31], [32]. In other words, we map 
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the SIR dynamics onto the bond percolation with open bond probabihty T. Strictly 
speaking, this mapping is not generally exact [331 IMl [35] . The transition point above 
which a global outbreak occurs and the mean final size of global outbreaks at a given 
T obtained from the analysis of the bond percolation are precise for the SIR model. 
However, the bond percolation fails to predict the probability of global outbreaks in the 
SIR model [33l [311 [35] . The correspondence between the bond percolation and the SIR 
model is exact if the infectious period in the SIR model is assumed to be constant. When 
analyzing the effects of defense strategies (Sec. 13. 2| ). we denote by Ti the probability that 
susceptible node i is infected by an infected neighbor. We set = T for all i unless 
otherwise stated. 

After a propagating attack, each node is either susceptible or removed (figure [T](a)). 
We call the network composed of removed nodes and that composed of susceptible nodes 
infected network and residual network, respectively. The infected network consists of 
a single component because an attack starts from a single seed. When T is below a 
first transition point Td, the size of the infected network is o{N), whereas the residual 
network includes a giant component of order 0{N). When < T < Tc2, where is a 
second transition point, both the infected and residual networks have a giant component 
of order 0{N). When T > Tc2, the giant component of the residual network is absent. 

We consider two defense strategies. In the random defense, Nf randomly selected 
nodes are vaccinated. In the degree-based defense, Nf nodes with the largest degrees 
are vaccinated. We assume that Tj decreases to Tj = Tf^ (0 < /v < 1) for a vaccinated 
node i (figure [D^b)). In epidemiology, the vaccine with /v = and that with /v 7^ are 
called perfect and leaky vaccination, respectively [221 [36]. 

3. Analysis 

3.1. SIR Model on Networks 

In the following, we consider uncorrelated, locally tree-like networks with arbitrary 
degree distributions. In this section, we quickly review the derivation of Td and in 
the absence of a defense strategy [211 [31] • 

We denote by pk the fraction of nodes with degree k. We denote by the 
distribution of the excess degree defined as the degree of the node reached by following 
a randomly selected link minus 1; qk = {k + l)pk+i/ (k), where (■) represents the average 
of a quantity weighted by pk- We define 




(1) 



k 




(2) 



The giant component of the infected network is present if [21] 
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Figure 2. Schematics of (a) F^{x), (b) G^{x,y), (c) P^ix,y), and (d) P^ix,y). 



To derive Tc2, let m be the probability that a node is not infected by a specific 
neighbor in a propagating attack. After a propagating attack, a node is susceptible 
with probability YlkPkU^ ~ ^o{u), and a node reached by following a link is susceptible 
with probability ^j^qkU^ = Fi{u). The following recursive relationship determines u: 

1-T + TFi(m). (4) 



u 



The probability that a node is susceptible in the residual network and its randomly 
selected neighbor is removed is equal to (1 — T)[l — Fi{u)] = u — Fi{u). Therefore, the 
probability that a node in the residual network has m susceptible neighbors, i.e., the 
degree distribution of the residual network, is represented as 



k=m 



Fi(nr[n-Fi(n)] 



k—m 



(5) 



The corresponding generating function Po{x) is given by 

^ OO OO / , 

TVl 2^ 2^ 
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1 r^/ N/ N Fr)\Fi(u)(x — 1) + u] , , 

Similarly, the generating function Pi{x) for the excess degree distribution of the residual 
network is represented as 

_ F,[F,{u){x - 1) + u] 

From Eq. ([7]), the mean excess degree of the residual network, denoted by (fccx)ros, is 
calculated as 

{ko.)rcs=^Pl{x) =F[{u). (8) 

ax x=i 

The giant component of the residual network is absent if T > Tc2, where Tc2 is 
determined by (fcex)res = F[iu) = 1 [2T| . 

5.^ Defense strategies against propagating attacks 

In the following, we calculate Td and Tc2 in the presence of the random or degree-based 
defense. 

3.2.1. Degree Distributions of Vaccinated and Unvaccinated Nodes in the Original 
Network Under a defense strategy, each node is either unvaccinated (U) or vaccinated 
(V) (figure [Hb)). First, we derive the degree distribution of U nodes and that of V 
nodes in the original network. Under both defense strategies, a randomly selected node 
is U with probability 1 — / and V otherwise. We denote the probability that a node 
reached by following a randomly selected link is a V node by /. Then, the probability 
that a node reached by following a randomly selected link is a U node is equal to 1 — /. 
Under the random defense, we obtain 

/ = /• (9) 
Under the degree-based defense, we obtain 

oo 

where /Ccut is the minimum degree of V nodes determined by 

oo 

Let and p^ be the degree distributions of U nodes and V nodes, respectively, 
and and be the distributions of the excess degrees of U nodes and V nodes, 
respectively. For the random defense, we obtain 

{p]^,pli} = {Pk,Pk}, (12) 
Uk^Qk] = Uk^Qk}, (13) 
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for all k. For the degree-based defense, we obtain 

r^u V. _ / {Pk/il - /),0} (A; < ^cut), .^ 

^^^'^'^^"l {0,W/} ihu.<k), ^''^ 

,u vi _ / {9fc/(l - /).0} {k < kcut), 



3.2.2. Derivation of Td We analyze the interaction of the generating functions for U 
nodes and those for V nodes. Related approaches have been adopted to analyze the 
percolation in networks that contain multiple types of nodes [371 [38] . We define the 
following generating functions 

oo 

F^{x)^Y.P'>^x^ (16) 

A;=0 
oo 

F^{x)^Y.P^,x', (17) 

fc=0 
oo 

Fnx)^J2l^x', (18) 

k=0 
oo 

Fnx)^J2^Ix'. (19) 

fc=0 

A schematic of Fq{x) is shown in figure [2]^a). 

We extend the theory developed in [39] to derive Td for our model. Consider the 
mean number Ni of removed nodes at distance i from the seed node. We obtain 

oo k / J \ 

^1 = T.P'^ J2[Jl^T + {k- m)TA](l - fTf'-'^ 

k=0 m=0 ^ ^ 



(t— + T/.— )Fo[(1 - f)x + fy]l^^^ = [T{1 - /) + TfJ]F^{l), (20) 



/ 1 \ 



k 

m f'k—m 



fTf 



k=0 m=0 

oo 



m' rk'—m' 



X {mT J2 Qk' E (10 f""'^ + ~ ^')Tm - fr'f 

k'=0 m'=0 ^ ^ 

oo ^' / h' \ 

+ {k-m)Tf^Y.'l^' E (^J[^'T + (fc'-m')TM(l-/r' 

A;'=0 m'=0 ^ 

F^(1)[T(1 - /) +T/J][r(l - f)F^\l)+TfjF^\l)l (21) 



and 

iv, = F^(i)[r(i-/) + r/./][T(i-/)FU'(i) + r/./F7'(i)ri= (^) iVi, (22) 
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for any i > 1. The infected network is a giant component if and only if N2 > Ni, or 
equivalently, 

^ = T(l - />u'(i) + TfJ~F^\l) > 1. (23) 

Therefore, T^i is determined from 

r,i[(l - />u'(i) + fjF^'il)] = 1. (24) 

In the case of the random defense, where (1) = (1) = F[{1) and f = f , 
Eq. (!24|) is reduced to 

Tci[(l-/) + /v/]^^^^^ = l. (25) 

3.2.3. Derivation ofTc2 The probabihty that a randomly selected U node has exactly 
m U neighbors is equal to Xlfclm^fc (m) (■'- ~ /)'"/'^~'". The generating function Gq (x, y) 
for this probability distribution, which is schematically shown in figure HJ^b), is given by 

GU(x,y) = J^pUj^ f j(l-/)-/^— xV""^ = FoU[(l-/>+/y].(26) 

The generating function G^{x^y) for the probability distribution that a node has m 
U neighbors is given by 

k=0 m=0 ^ 



Similarly, we obtain 



00 k / 7 \ 

G'^ix^y) = J^g^ E (){l-fr'f'-"'x^y'-'^ = FU[(l-/>+/y], (28) 

k=0 m=0 ^ ^ 

GUx,y) ^ ^gV^ j(l-/)-/'=-a;V- = F7[(l-/>+/y].(29) 



k=0 m=0 



Let MxY (-^, ^ € {^5 ^}) be the probability that an X node is not infected by 
a specific Y neighbor in a propagating attack. Then, a U node is susceptible with 
probability Gq (uuu^^uv), and a V node is susceptible with probability ^^(uvu^^vv)- 
A U node connected to a randomly selected link is susceptible with probability 
, u\jy) , and a V node connected to a random selected link is susceptible with 
probability CYiuyu, uyy). Then, u\ju, muv, uyu, and Myy satisfy the following recursive 
relationships: 

Muu = 1 - T + TG^{uuu, Muv), (30) 

Ujjy = 1 -T + TG'^{uYu,UYy), (31) 

Mvu = 1 - Tf^ + Tf^G^{uuu, Muv), (32) 

MVV = 1 - T/v + Tf^G'^{uyv, Mvv). (33) 



Robustness of networks against propagating attacks under vaccination strategies 9 

To derive the degree distribution of U nodes in the residual network, consider a 
focal U node. When the focal U node is adjacent to a U node, the probability that the 
U neighbor is (infected and) eventually removed and the U focal node is not infected by 
this U neighbor in a propagating attack is equal to 

(1 - T)[l - G^{uuu,uuy)] = uuu - G^{uuv,uuy). (34) 

When the focal U node is adjacent to a V node, the probability that the V neighbor is 
eventually removed and the focal U node is not infected by this V neighbor is equal to 

(1 - T)[l - G'({uYu,UYy)] = nuv - Gi{uyu,UyY). (35) 

Then, the probability that a U node has m susceptible U neighbors and n susceptible 
V neighbors and is not infected is equal to 

k=m+nm'=m ^ ^ 

^/k- m'\ ^v^^^^^ nvv)'^Kv - ^^(nvu, nw)]'""^'-". 



n 



Therefore, the generating function Pq (x, y) for the degree distribution of U nodes in 
the residual network (figure [2t^c)) is given by 

k 



\m' rk—m' 



^\m) muv)"'[muu - G^{uuu, muv)] 

m=0 ^ ^ 
k-m' _ m'\ 

^ ( jG'y(Mvu,Mvv)"[Muv - G'7(mvu,mvv)] 

n—n \ ^ J 
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We can similarly derive the degree distribution of V nodes in the residual network. 
When a focal V node is adjacent to a U node, the probability that the U neighbor is 
eventually removed and the focal V node is not infected by this U neighbor is equal to 

(1 - T/v)[l - G[^(muu,muv)] = Mvu - G'1^(muu,muv)- (37) 

When the focal V node is adjacent to a V node, the probability that the V neighbor is 
eventually removed and the focal V node is not infected by this V neighbor is equal to 

(1 - T/^)[l - GI{uyv,uyy)] = uvY - Gy(uvu,^ivv)- (38) 

The generating function PQ{x,y) for the degree distribution of V nodes in the residual 
network (figure EJ^d)) is given by 

Poi^. y) = rvr^ — -.f^pl E { ^) (1 - fr'I'-' 

G^(«vu,nvv)^ ^,\m'J 
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X )g1^Ku,wuv)"^Ku - G^Ku, wuv)]"''-"'^'" 

m=0 ^ ^ 

Gy(uvu,«vv)"Kv - Gy(nvu,«vv)]'-"^'-"y" 

»i— n \ ^ J 

(39) 



n=0 

^^[^^(uuu, ^iuv)(a; - 1) + ^ivu, ^^(mvu, UYv){y - 1) + ^^vv] 



G'^(uvu,'Uvv) 

Similarly, the generating functions for the excess degree distributions of U nodes 
and V nodes in the residual network are given by 
pU/ N _ GY[GY{u\ju,uuy){x - 1) + ^uu,gy(nvu,nvv)(y - 1) + ^uv] 

f-^i («uu,«uv) 

and 

^ _ Gy[gy(Muu,Muv)(x - 



pV(^^ ^) ^ GX[GY{u\ju,uuy){x - 1) +Mvu,Gy(nvu,«vv)(y - 1) + Uyv] ^^^-^ 



respectively. 

When we follow a randomly selected link in the residual network, we reach a U 
node with probability 

(1 - /)Gy(MuU,Muv) 

(1 - /)Gy(MuU,Muv) + /Gy(Mvu,Mvv) 

Otherwise, we reach a V node. Then, the mean excess degree of the residual network is 
given by 

^^^'•'^^ \dx dy) (1 - f)GY{uu^, «uv) + /Gr(«vu, «vv) 

= (1 - f)F^'[{l - f)u^^ + /nuv] + />r[(l - /)nvu + fUyy]. (42) 
Because the second transition point is given by (/Cex)res = 1? we find Tc2 by solving 

(1 - /)FU'[(1 - />uu + /nuv] + />i'''[(l - f)uyv + hw] = 1. (43) 

3.2.4- Component Size The generating function formalism also gives the component 
size of the infected and residual networks. Because a randomly selected U node is 
susceptible with probability Gq (uuu, ^^uv) and a randomly selected V node is susceptible 
with probability Gq (mvu; ''^vv); the component size of the infected network Si^{ is 
represented as 

5i„f = 1 - (1 - /)GU(nuu, nuv) - fG^iuyv, ^w). (44) 

To derive the largest component size of the residual network, we denote by Hq{x) 
and Hq{x) the generating functions for the size of the finite components of the residual 
networks to which a randomly selected U node and V node belong [3] . The average size 
Srcs of the largest component of the residual network relative to the size of the original 
network is given by 

G'U(nuu,«uv)i^o"(l) + G'^Ku,nvv)ifo''(l) 



'S'res — (1 'S'inf) 



G^(muu, muv) + G^(mvu, Mvv) 



(45) 



Robustness of networks against propagating attacks under vaccination strategies 



11 



where 

H^{1) = P:^{AB), (46) 
H^{l) = P^{A,B), (47) 

A^H]^{l)=P^iA,B), (48) 
B = H^{1) = P^{A,B). (49) 

3.2.5. Degree- dependent Probability of Infection We denote by U{k) the probabihty 
that a node with degree k is uninfected. Because a U node is susceptible with probabihty 
Gq (uuu) ^uv) and a V node is susceptible with probability Gq^uyu, uyy), U{k) is given 
by 



Y,VkU{k) = (1 - /)GU(muu,muv) + /G^(wvu, 



MVVj 



fc=0 



= (1 - /) 5^|>I^[(1 - />uu + fu^jyf + / J]pr[(l - /)«vu + /«vv]'. (50) 

fc=0 fc=0 

For the random defense, p^ and p^ are given by Eq. (fl2l) . Therefore, Eq. fl50l) is reduced 
to 

oo oo 

5^PA:t/(A:) = ^R.{(1 - /)[(1 - /)muu + /wuv]' + /[(I - /)wvu + fu^vf}. (51) 

fc=0 fc=0 

Because Eq. ( |5T1) is satisfied for an arbitrary degree distribution {pk}, we obtain 

?7(A:) = (1 - /)[(1 - /)wuu + /wuv]' + /[(I - /)wvu + fuy^f. (52) 

For the degree-based defense, p^ and are given by Eq. ^\M . Therefore, Eq. fl50|l is 
reduced to 

oo fccut-1 OO 

^PkU{k) = pfc[(l - />uu + fuvvf + Yl ^'^[(^ ~ ^"^''vu + fuYvf. (53) 

fc=0 k=0 k=kcut 

Because Eq. fl53l) is satisfied for arbitrary {pfc}, we obtain 

= / [(^ " /)^uu + /^uv]'' (A: < /^cut), 

^ ^ ^ [(1 _ + /«vv]' (fccut<A;). ^ ^ 



4. Numerical Results 

On the basis of the theoretical results obtained in Sec. [3], we numerically examine the 
robustnesses of different networks against the propagating attack. We consider the 
following three combinations of networks and defense strategies: 

(i) random graph with random defense, 

(ii) uncorrelated SF network with random defense, and 

(iii) uncorrelated SF network with degree-based defense. 
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Figure 3. Component size 5i„f of the infected network (dashed hnes) and the largest 
component size S^cs of the residual network (solid lines) for (a) random graph with the 
random defense (case (i)), (b) SF network with the random defense (case (ii)), and (c) 
SF network with the degree-based defense (case (iii)). We set / = 0.05. The different 
lines correspond to /v = 0.00, 0.50, and 1.00, from right to left. 



We assume that the SF network has degree distribution oc k~^, where k > /cmin = 2, 
which yields (k) ^ 3.19. For a fair comparison, we assume that (k) of the random 
graph is equal to that of the SF network. For the random graph, we obtain Fq{x) = 
Fi{x) = exp[{k){x — 1)]- For the SF network, we obtain Fq{x) = Yll^=oPk^^ — 
Li3(x)/Li3(l), Fi{x) = Yl'k=o^PkX^~^ / Y.'k=o^Pk = Li2(x)/xLi2(l), where Li„(x) is the 
nth polylogarithm of x, i.e., Li„(a;) = YlT=i x^k~^. In the random defense (cases (i) and 
(ii)), Fq{x) = Fq{x) = Fo{x) and F^{x) = F^{x) = Fi{x). In the degree-based defense 
(case (iii)), we vaccinate all the nodes whose degree is equal to or larger than /ccut such 
that the fraction of V nodes does not exceed /. Then, we vaccinate some nodes with 
degree kcut — 1 such that the total fraction of V nodes is equal to /. The other nodes 
with degree fccut — 1 are U nodes. The generating functions for U and V nodes are given 
by 

^oi^) = Y.P^'''' ~ (55) 



fc=0 

k' 



F^{x) = Y,Qkx'-Cqk'^y^^\ (56) 

fc=0 
oo 

OO 

F7{^)= J2 ^kX^ + Cqy^y^^\ (58) 

^ — fccut 

where C = {f — Yl'kLk^^t'P^'^ I^Kut ^cut = ^cut — 1- For each case, we apply the 
results of Sec. 13.21 to obtain the component sizes and transition points. For given / and 
/v, we substitute F^{x), F^{x), F^{x), and F^{x) in Eq. to obtain T^. We also 
substitute Fq{x), F^{x), Fq{x), and F^{x) in Eqs. fl30|) - fl33|l to determine u\j\j, muv, 
uyij, and uyy. Then, we use the obtained values of uuu, muv, mvuj and uyy to derive 
Tc2, 5'inf, and S^ef from Eqs. (143|1 . (l44l) . and (l45ll . respectively. 
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Figure 4. Results for the random graph with the random defense (case (i)). (a) Td 
and (b) Tc2 as functions of / and /v The solid lines indicate Td = 1 in (a) and 
Tea = 1 in (b). The dashed lines indicate /(I - /v) = 0.1 and /(I - f^) = 0.2. (c) Tc2 
conditioned by /(I — fy) = 0.1 (circles) and 0.2 (diamonds). 



We first examine the component sizes ^inf of tlie infected network and S^cs of the 
residual network using Eqs. (jS]) and f HSjl . The results for case (i) with / = 0.05 are 
shown in figure |3]^a). When T is small, S'inf = (in the limit of the infinite network 
size) holds true irrespective of the value of /v, as shown by different dashed lines in 
figure Et^a). When T > Td ~ 0.33, S'inf increases and S^cs decreases with an increase in 
T. When T > ~ 0.6, we obtain Sj-cs = 0. The precise values of Td and Tc2 vary with 
/v. Efficient vaccines (i.e., small /v) make S'inf small and S'rcs large and increase Td and 

Te2. 

The results for case (ii) are shown in figure [Ht^b). Td is much smaller than in case 
(i); Eq. ( j25l) implies Td = when FY'{1) diverges. For /v 7^ 0, Srcs behaves qualitatively 
the same as that for case (i). In particular, the value of Tc2 for = 1 is only slightly 
smaller than that for case (i). For = 0, S'rcs behaves differently; S^es > even when 
T = 1, implying that the residual network is never destroyed. 

The results for case (iii) are shown in figure EJj^c). We find that the degree- 
based defense is effective at decreasing (increasing) Sinf (Sj-cs) as compared with the 
random defense (case (ii)). In particular, the perfect vaccines (i.e., /v = 0) contain the 
propagating attack when T < Td ~ 0.63. We find from Eq. ( l24l) that Td 7^ only 
for /v = 0. Therefore, any leaky vaccination (i.e., /v > 0) fails to contain the attack 
because F^'{1) diverges. We remark that, for /v = 0, Eq. ( IMl) is reduced to a known 
result [26] 

We also remark that, for = 1, the value of Tc2 coincides with that for case (ii) as 
expected; the effect of vaccination is completely null when fy = 1. 

In figures HI and [HI Td and Tc2 obtained from Eqs. fl24p and fH51) are shown for 
the three cases. 
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Figure 5. Results for the SF network with random defense (case (ii)). (a) Tc2 as 
a function of / and f^. Tlie solid line indicates = 1- The dashed lines indicate 
/(I - /^) = 0.1 and /(I - /v) = 0.2. (b) T^a conditioned by /(I - /v) = 0.1 (circles) 
and 0.2 (diamonds). 
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Figure 6. Results for the SF network with the degree-based defense (case (iii)). (a) 
Tc2 as a function of / and fv The solid line indicates Tc2 = 1- The dashed lines 
indicate /(I - ,U) = 0.1 and /(I - /v) = 0.2. (b) Tc2 conditioned by /(I - /v) = 0.1 
(circles) and 0.2 (diamonds). 



Tci monotonically increases with / and monotonically decreases with in case (i) 
(figure Hl^a)). When Td = 1 (the red region bounded by the bold fine in figure Hl^a)), 
a global infection does not occur at any infection rate. On the basis of Eq. ( IMl) . the 
boundary of the region T^i = 1 is given by 

In case (ii), Td = for any / and unless (/, /v) = (1,0). In case (iii), = if 
/v > or (/, /v) = (0, 0) for the following reason. If > 0, ^^'(l) in Eq. ([23]) diverges. 
If / = 0, FU'(1) in Eq. ([23]) diverges (and ^^'(1) = 0). 

Similar to Td, Tc2 monotonically increases with / and monotonically decreases with 
/v in all the cases (figures ID^b), |5]^a), and[6]^a)). When Tc2 = 1, the residual network 
is not disintegrated at any infection rate (the red regions bounded by the bold line in 
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Figure 7. Probability U{k) tliat a node witli degree k in the SF network is uninfected, 
(a) Random defense (case (ii)). (b) Degree-based defense (case (iii)). We set To = 0.4 
and / = 0.05. The daslied line in (b) indicates fccut = 8. 



figures |4](b), [5]^a), and[6]^a)). 

When /v = 0, the SF network yields Tc2 = 1 for any / > 0, either for the random 
defense (figure [5]|^a)) and the degree-based defense (figure [6](a)). In this situation, V 
nodes are never infected such that they percolate in a SF network. This phenomenon is 
identical to the standard site percolation on SF networks with 7 < 3; randomly located 
nodes percolate even at an infinitesimally small occupation probability flU\ [TTj. 

The comparison of the values of Tc2 among the three cases based on figures Hl^b), 
[5]^a), and[6](a) suggests the following. When the vaccines are not efficient (i.e., ~ 1), 
the network is the most robust for case (i) among the three cases in the sense that Tc2 
is the largest. In this situation, Tc2 for the SF network is slightly smaller than that for 
the random graph (also see figure [3]). When takes an intermediate value, the network 
is the most robust for case (iii). The network is the most fragile for case (ii) except for 
fairly small /v and /, in which case the residual network mainly composed of vaccinated 
nodes can percolate on the SF network and not on the random graph. 

In reality, there may be a trade-off between the required number of vaccines and the 
efficiency of vaccines. Therefore, we measure Tc2 under the constraints /(I — /v) = 0.1 
and /(I — /v) = 0.2; / and 1 — /v represent the number and efficiency of vaccines, 
respectively. As shown in figures IH^c),[5](b), and[6](b), the network becomes progressively 
robust as / decreases (i.e., small number of vaccines) and /v decreases (i.e., high 
efficiency of vaccines) in all the three cases. 

The difference between the random and degree-based defenses in the SF network is 
apparent when we look at U{k), i.e., the probability that a node with degree k remains 
uninfected. In figures [Tl^a) and[7](b), U{k) values are shown for cases (ii) (Eq. fl52l) ) 
and (iii) (Eq. ( IMI) ). respectively. As expected, the degree-based defense cuts down the 
probability of infection for nodes with high degrees, whereas the random defense does 
not. Furthermore, the indirect effect of vaccines, i.e., the suppression of the infection 
probability for U nodes, is strong under the degree-based defense. For example, some 
fractions of U nodes with degree k < /Ccut (below the dashed line in figure [7] (b)) remain 
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uninfected. 



5. Summary 

We investigated the effect of defense strategies on the robustness of networks against 
propagating attacks by extending the generating function framework developed by 
Newman [21]. In particular, under vaccination strategies, we analytically obtained the 
second critical infection rate above which the residual network is disintegrated. We 
applied the analytical results to the three cases: the random graph with the random 
defense, the SF network with the random defense, and the SF network with the degree- 
based defense. The first critical point of the SF network depends crucially on whether 
the effect of the vaccines is perfect or leaky. Even under the degree-based defense, the 
leaky vaccines cannot prevent a global outbreak on heterogeneous networks. We also 
found that, under a trade-off between the number and efficiency of vaccines, it is better 
to administer a small number of vaccines with high efficiency than vice versa to prevent 
the residual network from being disintegrated. 
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